RECENT RESULTS 
ON THE PERIODIC LORENTZ GAS 



FRANgOIS GOLSE 

Abstract. The Drude-Lorentz model for the motion of electrons 
in a solid is a classical model in statistical mechanics, where elec- 
trons are represented as point particles bouncing on a fixed system 
of obstacles (the atoms in the solid) . Under some appropriate scal- 
ing assumption — known as the Boltzmann-Grad scaling by anal- 
ogy with the kinetic theory of rarefied gases — this system can be 
described in some limit by a linear Boltzmann equation, assuming 
that the configuration of obstacles is random [G. Gallavotti, [Phys. 
Rev. (2) 185 (1969), 308]). The case of a periodic configuration 
of obstacles (like atoms in a crystal) leads to a completely differ- 
ent limiting dynamics. These lecture notes review several results 
on this problem obtained in the past decade as joint work with J. 
Bourgain, E. Caglioti and B. Wennberg. 
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Introduction: from particle dynamics to kinetic models 

The kinetic theory of gases was proposed by J. Clerk Maxwell [3HI35] 
and L. Boltzmann [5] in the second half of the XlXth century. Because 
the existence of atoms, on which kinetic theory rested, remained con- 
troversial for some time, it was not until many years later, in the XXth 
century, that the tools of kinetic theory became of common use in var- 
ious branches of physics such as neutron transport, radiative transfer, 
plasma and semiconductor physics... 

Besides, the arguments which Maxwell and Boltzmann used in writ- 
ing what is now known as the "Boltzmann collision integral" were far 
from rigorous — at least from the mathematical viewpoint. As a matter 
of fact, the Boltzmann equation itself was studied by some of the most 
distinguished mathematicians of the XXth century — such as Hilbert 
and Carleman — before there were any serious attempt at deriving 
this equation from first principles (i.e. molecular dynamics.) Whether 
the Boltzmann equation itself was viewed as a fundamental equation 
of gas dynamics, or as some approximate equation valid in some well 
identified limit is not very clear in the first works on the subject - 
including Maxwell's and Boltzmann's. 

It seems that the first systematic discussion of the validity of the 
Boltzmann equation viewed as some limit of molecular dynamics - 
i.e. the free motion of a large number of small balls subject to binary, 
short range interaction, for instance elastic collisions — goes back to 
the work of H. Grad [26]. In 1975, O.E. Lanford gave the first rigorous 
derivation of the Boltzmann equation from molecular dynamics - 
his result proved the validity of the Boltzmann equation for a very short 
time of the order of a fraction of the reciprocal collision frequency. (One 
should also mention an earlier, "formal derivation" by C. Cercignani 
[12] of the Boltzmann equation for a hard sphere gas, which consider- 
ably clarified the mathematical formulation of the problem.) Shortly 
after Lanford's derivation of the Boltzmann equation, R. Illner and M. 
Pulvirenti managed to extend the validity of his result for all positive 
times, for initial data corresponding with a very rarefied cloud of gas 
molecules [2"T] . 
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An important assumption made in Boltzmann's attempt at justifying 
the equation bearing his name is the "Stosszahlansatz" , to the effect 
that particle pairs just about to collide are uncorrelated. Lanford's ar- 
gument indirectly established the validity of Boltzmann's assumption, 
at least on very short time intervals. 

In applications of kinetic theory other than rarefied gas dynamics, 
one may face the situation where the analogue of the Boltzmann equa- 
tion for monatomic gases is linear, instead of quadratic. The linear 
Boltzmann equation is encountered for instance in neutron transport, 
or in some models in radiative transfer. It usually describes a situa- 
tion where particles interact with some background medium — such 
as neutrons with the atoms of some fissile material, or photons subject 
to scattering processes (Rayleigh or Thomson scattering) in a gas or a 
plasma. 

In some situations leading to a linear Boltzmann equation, one has 
to think of two families of particles: the moving particles whose phase 
space density satisfies the linear Boltzmann equation, and the back- 
ground medium that can be viewed as a family of fixed particles of 
a different type. For instance, one can think of the moving particles 
as being light particles, whereas the fixed particles can be viewed as 
infinitely heavier, and therefore unaffected by elastic collisions with the 
light particles. Before Lanford's fundamental paper, an important - 
unfortunately unpublished — preprint by G. Gallavotti [T5] provided 
a rigorous derivation of the linear Boltzmann equation assuming that 
the background medium consists of fixed, independent like hard spheres 
whose centers are distributed in the Euclidian space under Poisson's 
law. Gallavotti's argument already possessed some of the most remark- 
able features in Lanford's proof, and therefore must be regarded as an 
essential step in the understanding of kinetic theory. 

However, Boltzmann's Stosszahlansatz becomes questionable in this 
kind of situation involving light and heavy particles, as potential cor- 
relations among heavy particles may influence the light particle dy- 
namics. Gallavotti's assumption of a background medium consisting of 
independent hard spheres excluded this this possibility. Yet, strongly 
correlated background media are equally natural, and should also be 
considered. 

The periodic Lorentz gas discussed in these notes is one example of 
this type of situation. Assuming that heavy particles are located at the 
vertices of some lattice in the Euclidian space clearly introduces about 
the maximum amount of correlation between these heavy particles. 
This periodicity assumption entails a dramatic change in the structure 
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FIGURE 1. Left: Paul Drude (1863-1906); right: Hen- 
drik Antoon Lorentz (1853-1928) 



of the equation that one obtains under the same scaling limit that 
would otherwise lead to a linear Boltzmann equation. 

Therefore, studying the periodic Lorentz gas can be viewed as one 
way of testing the limits of the classical concepts of the kinetic theory 
of gases. 

Acknowledgements. 

Most of the material presented in these lectures is the result of col- 
laboration with several authors: J. Bourgain, E. Caglioti, H.S. Dumas, 
L. Dumas and B. Wennberg, whom I wish to thank for sharing my 
interest for this problem. I am also grateful to C. Boldighrini and G. 
Gallavotti for illuminating discussions on this subject. 

1. The Lorentz kinetic theory for electrons 

In the early 1900's, P. Drude [16] and H. Lorentz [30] independently 
proposed to describe the motion of electrons in metals by the methods 
of kinetic theory. One should keep in mind that the kinetic theory of 
gases was by then a relatively new subject: the Boltzmann equation for 
monatomic gases appeared for the first time in the papers of J. Clerk 
Maxwell |35j and L. Boltzmann [5]. Likewise, the existence of electrons 
had been established shortly before, in 1897 by J.J. Thomson. 
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The basic assumptions made by H. Lorentz in his paper [30J can be 
summarized as follows. 

First, the population of electrons is thought of as a gas of point 
particles described by its phase-space density / = f(t,x,v), that is the 
density of electrons at the position x with velocity v at time t. 

Electron-electron collisions are neglected in the physical regime con- 
sidered in the Lorentz kinetic model — on the contrary, in the classical 
kinetic theory of gases, collisions between molecules are important as 
they account for momentum and heat transfer. 

However, the Lorentz kinetic theory takes into account collisions 
between electrons and the surrounding metallic atoms. These collisions 
are viewed as simple, elastic hard sphere collisions. 

Since electron-electron collisions are neglected in the Lorentz model, 
the equation governing the electron phase-space density / is linear. 
This is at variance with the classical Boltzmann equation, which is 
quadratic because only binary collisions involving pairs of molecules 
are considered in the kinetic theory of gases. 

With the simple assumptions above, H. Lorentz arrived at the fol- 
lowing equation for the phase-space density of electrons / = f(t,x,v): 



In this equation, C is the Lorentz collision integral, which acts on 
the only variable v in the phase-space density /. In other words, for 
each continuous function <p = <f)(v), one has 



The other parameters involved in the Lorentz equation are the mass 
m of the electron, and N at , r at respectively the density and radius of 
metallic atoms. The vector field F = F(t, x) is the electric force. In the 
Lorentz model, the self-consistent electric force — i.e. the electric force 
created by the electrons themselves — is neglected, so that F take into 
account the only effect of an applied electric field (if any). Roughly 
speaking, the self consistent electric field is linear in /, so that its con- 
tribution to the term F • V«/ would be quadratic in /, as would be any 
collision integral accounting for electron-electron collisions. Therefore, 
neglecting electron-electron collisions and the self-consistent electric 
field are both in accordance with assuming that /<1. 



(d t + v-V x + ±F{t, x) ■ V«)/(t, x, v) = N at r 2 at \v\C{f)(t, x, v) . 




and the notation 



C(f)(t, x, v) designates C(f(t, x, -))(v) . 
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Figure 2. The Lorentz particle path 

The line of reasoning used by H. Lorentz to arrive at the kinetic 
equations above is based on the postulate that the motion of electrons 
in a metal can be adequately represented by a simple mechanical model 
- a collisionless gas of point particles bouncing on a system of fixed, 
large spherical obstacles that represent the metallic atoms. Even with 
the considerable simplification in this model, the argument sketched in 
the article [SU] is little more than a formal analogy with Boltzmann's 
derivation of the equation now bearing his name. 

This suggests the mathematical problem, of deriving the Lorentz ki- 
netic equation from a microscopic, purely mechanical particle model. 
Thus, we consider a gas of point particles (the electrons) moving in 
a system of fixed spherical obstacles (the metallic atoms). We as- 
sume that collisions between the electrons and the metallic atoms are 
perfectly elastic, so that, upon colliding with an obstacle, each point 
particle is specularly reflected on the surface of that obstacle. 

Undoubtedly, the most interesting part of the Lorentz kinetic equa- 
tion is the collision integral which does not seem to involve F. Therefore 
we henceforth assume for the sake of simplicity that there is no applied 
electric field, so that 

F(t,x) = 0. 

In that case, electrons are not accelerated between successive col- 
lisions with the metallic atoms, so that the microscopic model to be 
considered is a simple, dispersing billiard system — also called a Sinai 
billiard. In that model, electrons are point particles moving at a con- 
stant speed along rectilinear trajectories in a system of fixed spherical 
obstacles, and specularly reflected at the surface of the obstacles. 

More than 100 years have elapsed since this simple mechanical model 
was proposed by P. Drude and H. Lorentz, and today we know that 
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the motion of electrons in a metal is a much more complicated physical 
phenomenon whose description involves quantum effects. 

Yet the Lorentz gas is an important object of study in nonequilibrium 
satistical mechanics, and there is a very significant amount of literature 
on that topic — see for instance [H] and the references therein. 

The first rigorous derivation of the Lorentz kinetic equation is due 
to G. Gallavotti [TBI tlH], wri ° derived it from from a billiard system 
consisting of randomly (Poisson) distributed obstacles, possibly over- 
lapping, considered in some scaling limit — the Boltzmann-Grad limit, 
whose definition will be given (and discussed) below. Slightly more 
general, random distributions of obstacles were later considered by H. 
Spohn in |^3] . 

While Gallavotti's theorem bears on the convergence of the mean 
electron density (averaging over obstacle configurations), C. Boldrigh- 
ini, L. Bunimovich and Ya. Sinai |4] later succeeded in proving the 
almost sure convergence (i.e. for a.e. obstacle configuration) of the 
electron density to the solution of the Lorentz kinetic equation. 

In any case, none of the results above says anything on the case of a 
periodic distribution of obstacles. As we shall see, the periodic case is 
of a completely different nature — and leads to a very different limiting 
equation, involving a phase-space different from the one considered by 
H. Lorentz — i.e. R 2 xS 1 — on which the Lorentz kinetic equation is 
posed. 

The periodic Lorentz gas is at the origin of many challenging math- 
ematical problems. For instance, in the late 1970s, L. Bunimovich and 
Ya. Sinai studied the periodic Lorentz gas in a scaling limit different 
from the Boltzmann-Grad limit studied in the present paper. In [7], 
they showed that the classical Brownian motion is the limiting dynam- 
ics of the Lorentz gas under that scaling assumption — their work was 
later extended with N. Chernov: see [8]. This result is indeed a major 
achievement in nonequilibrium statistical mechanics, as it provides an 
example of an irreversible dynamics (the heat equation associated with 
the classical Brownian motion) that is derived from a reversible one 
(the Lorentz gas dynamics). 



2. The Lorentz gas in the Boltzmann-Grad limit 
with a Poisson distribution of obstacles 

Before discussing the Boltzmann-Grad limit of the periodic Lorentz 
gas, we first give a brief description of Gallavotti's result [TEIEI] for the 
case of a Poisson distribution of independent, and therefore possibly 



8 



F. GOLSE 



overlapping obstacles. As we shall see, Gallavotti's argument is in some 
sense fairly elementary, and yet brilliant. 

First we define the notion of a Poisson distribution of obstacles. 
Henceforth, for the sake of simplicity, we assume a 2-dimensional set- 
ting. 

The obstacles (metallic atoms) are disks of radius r in the Euclidian 
plane R 2 , centered at ci, c 2 , . . . , Cj, . . . G R 2 . Henceforth, we denote by 

{c} — {ci,C2, ■ ■ ■ , Cj, • • •} = a configuration of obstacle centers. 

We further assume that the configurations of obstacle centers {c} 
are distributed under Poisson's law with parameter n, meaning that 

Prob({{c} | #(A n {c}) = p}) = e-^lM^! , 

where \A\ denotes the surface, i.e. the 2-dimensional Lebesgue measure 
of a measurable subset A of the Euclidian plane R 2 . 

This prescription defines a probability on countable subsets of the 
Euclidian plane R 2 . 

Obstacles may overlap: in other words, configurations {c} such that 

for some j ^ k G {1, 2, . . .}, one has |q — Cj\ < 2r 

are not excluded. Indeed, excluding overlapping obstacles means re- 
jecting obstacles configurations {c} such that |q — Cj\ < 2r for some 
i,j G N. In other words, Prob(ci{c}) is replaced with 

\ II l|ci- Cj -|>2rProb(d{c}) , 

i>j>0 

(where Z > is a normalizing coefficient.) Since the term 
1 1| l\ Cl - Cj \>2r is not of the form ]j0 fe (c fc ) , 

i>j>0 k>0 

the obstacles are no longer independent under this new probability 
measure. 

Next we define the billiard flow in a given obstacle configuration {c}. 
This definition is self-evident, and we give it for the sake of complete- 
ness, as well as in order to introduce the notation. 

Given a countable subset {c} of the Euclidian plane R 2 , the billiard 
flow in the system of obstacles defined by {c} is the family of mappings 

(X(t;;;{c}),V(t;;;{c})): (r 2 \ (Jb (c„ r)^ xS'O 

defined by the following prescription. 
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Whenever the position X of a particle lies outside the surface of any 
obstacle, that particle moves at unit speed along a rectilinear path: 

X(t; x, v, {c}) = V(t; x, v, {c}) , 

V(t; x, v, {c}) = , whenever \X(t; x, v, {c}) — q| > r for all i , 

and, in case of a collision with the i-th obstacle, is specularly reflected 
on the surface of that obstacle at the point of impingement, meaning 
that 

X(t + 0; x, v, {c}) = X(t - 0; x, v, {c}) G dB(a, r) , 

X(t; x, v, {c}) — c. 



V(t + 0;x,v, {c}) = U 



V(t-0;x,v,{c}), 



r 

where TZ[u] denotes the reflection with respect to the line (Rcj)- 1 : 

7?.[c<j]t> = v — 2(u ■ v)oj , |u;| = 1 . 

Then, given an initial probability density fV\ = fV£>(x,v) on the 
single-particle phase-space with support outside the system of obstacles 
defined by {c}, we define its evolution under the billiard flow by the 
formula 

/(*, x, v, {c}) = fl%(X(-t; x, v, {c}), V(-t; x, v, {c})) , t > . 

Let Ti(x, v, {c}), t 2 (x, v, {c}), . . . , Tj(x, v, {c}), ... be the sequence of 
collision times for a particle starting from x in the direction — v at t — 
in the configuration of obstacles {c}: in other words, 

T j( x ,v,{c}) = 

sup{t | #{s G [0,t] | dist(X(-s,x,v, {c}); {c}) = r} = j - 1} . 

Denoting r = and Ar fe = r k — r fc _ 1? the evolved single-particle 
density / is a.e. defined by the formula 

f(t, x, v, {c}) = f n (x - tv, v)l t<T1 

+ J2f n L-Jj\T k V(-T k )-(t - T^Vi-Tf), V(-T+) ) l Tj <t<r j+1 ■ 

j>i \ k=i J 

In the case of physically admissible initial data, there should be no 
particle located inside an obstacle. Hence we assumed that /?"■, = 
in the union of all the disks of radius r centered at the Cj G {c}. By 
construction, this condition is obviously preserved by the billiard flow, 
so that f(t, x, v, {c}) also vanishes whenever x belongs to a disk of 
radius r centered at any Cj G {c}. 
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As we shall see shortly, when dealing with bounded initial data, this 
constraint disappears in the (yet undefined) Boltzmann-Grad limit, as 
the volume fraction occupied by the obstacles vanishes in that limit. 

Therefore, we shall henceforth neglect this difficulty and proceed as 
if f m were any bounded probability density on R 2 x S 1 . 

Our goal is to average the summation above in the obstacle configu- 
ration {c} under the Poisson distribution, and to identify a scaling on 
the obstacle radius r and the parameter n of the Poisson distribution 
leading to a nontrivial limit. 

The parameter n has the following important physical interpretation. 
The expected number of obstacle centers to be found in any measurable 
subset VL of the Euclidian plane R 2 is 

^pProb({{c} | n {c}) = p}) = J2pe- nlnl ^p^- = n\Q\ 
so that 

n = # obstacles per unit surface in R 2 . 
The average of the first term in the summation defining f(t, x, v, {c}) 

is 

f m {x - tv, v)(l t<T1 ) = f m (x - tv, v)e' n2rt 

(where ( • ) denotes the mathematical expectation) since the condition 
t < T\ means that the tube of width 2r and length t contains no obstacle 
center. 

Henceforth, we seek a scaling limit corresponding to small obstacles, 
i.e. r — > and a large number of obstacles per unit volume, i.e. n — > oo. 
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There are obviously many possible scalings satisfying this require- 
ment. Among all these scalings, the Boltzmann-Grad scaling in space 
dimension 2 is defined by the requirement that the average over ob- 
stacle configurations of the first term in the series expansion for the 
particle density / has a nontrivial limit. 

Boltzmann-Grad scaling in space dimension 2 

In order for the average of the first term above to have a nontrivial 
limit, one must have 

r — > + and n — > +oo in such a way that 2nr — > o~ > . 



Under this assumption 

(f in (x - tv,v)l t<Tl ) - f in (x - tv,v)e-^ . 

Gallavotti's idea is that this first term corresponds with the solution 
at time t of the equation 

{pt + v • V x )f = —nrf / cos(v , uj)cLuj = —2nrf 



=i 

uj-v>0 



f £in 

\t=0 ~ •> 

that involves only the loss part in the Lorentz collision integral, and 
that the (average over obstacle configuration of the) subsequent terms 
in the sum defining the particle density / should converge to the 
Duhamel formula for the Lorentz kinetic equation. 

After this necessary preliminaries, we can state Gallavotti's theorem. 

Theorem 2.1 (Gallavotti [19]). Let f m be a continuous, bounded prob- 
ability density on R 2 x S 1 , and let 

f r {t, X, V, {C}) = f n ((X r , V r )(~t, X, V, { C })) , 

where (t,x,v) i— > (A r , V r )(t, x, v , {c}) is the billiard flow in the system 
of disks of radius r centered at the elements of {c}. Assuming that 
the obstacle centers are distributed under the Poisson law of parameter 
n = a/2r with a > 0, the expected single particle density 

(f r (t,x,v,-)) -> f(t,x,v) m L\R 2 x S 1 ) 

uniformly on compact t-sets, where f is the solution of the Lorentz 
kinetic equation 

(d t + v ■ V x )f + af = a fit, x, R[(3]v) sin 

Jo 

£ I _ fin 

J \t=0 ~ J ' 



2 4 ' 
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Figure 4. The tube T(t, c\,C2) corresponding with the 
third term in the series expansion giving the particle den- 
sity 



where R[/3] denotes the rotation of an angle (3. 

End of the proof of Gallavotti's theorem. The general term in the sum- 
mation giving f(t,x,v, {c}) is 

f in (^-i^r k V r (-r k )-(t - T j )V(-T+) 1 V r (-T+)\ l Tj <t<r j+1 , 
and its average under the Poisson distribution on {c} is 



I f m k-E &r k V r (-T k ) -(t- r y ) V(-r+), V 
J \ k=l 

r -n\T(t;cx,...,c,)\ njdc l--- dc 3 

r- 

where T{t; ci, . . . ,Cj) is the tube of width 2r around the particle tra- 
jectory colliding first with the obstacle centered at C\, . . . , and whose 
j-th collision is with the obstacle centered at cj. 
As before, the surface of that tube is 

\T(t; Cl ,..., Cj )\ =2rt + 0{r 2 ). 

In the j-th term, change variables by expressing the positions of the j 
encountered obstacles in terms of free flight times and deflection angles: 

(ci,...,^) i-> (n,...,^-; Pi,..., (3j). 

The volume element in the j-th integral is changed into 

**^=r>'siii& ...sin | ^ . . .^-dr x . . .d Tj . 
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Figure 5. The substitution (c 1 ,c 2 ) h-> (n, r 2 , ft, ft) 

The measure in the left-hand side is invariant by permutations of 
ci, . . . ,Cj] on the right-hand side, we assume that 

Ti < r 2 < . . . < Tj , 

which explains why factor disappears in the right-hand side. 

The substitution above is one-to-one only if the particle does not hit 
twice the same obstacle. Define therefore 
A r (T, x, v) = {{c} | there exists < t\ < t 2 < T and j'gN s.t. 

dist(X r (t!, x, v, {c}), Cj) = dist(X r (t 2 , x, v, {c}), c 3 ) = r} 

= [J{{c} | dist(X r (t,:r,i;, {c}),Cj) = r for some < ti < t 2 < T} , 



i>i 



and set 



f™(t, X, V, {c}) = f r (t, X, V, {c}) - f r H (t, X, V, {c}) , 



fr (t,X,V, {c}) = f r (t,X,V,{c})l Ar (T,x,v)({c}) , 

respectively the Markovian part and the recollision part in f r . 

After averaging over the obstacle configuration {c}, the contribution 
of the j-th term in f^ 1 is, to leading order in r: 



(2nr) J e 



j -2nrt 



I I 

J0<T1 <...<Ti<t J\0, 



sin 



Pi 



. sin 



ft d/3i 



xrlx-J2^r k R 

V k=l 

It is dominated by 



0<T 1 <...<T j <t J[0,2np 
k-l 

5> 



2 4 



dri . . . drj 



i=i 



V-{t-Tj)R 



3-1 

5> 



z=i 



j=i 



0((7) J e- 0(CT)t 
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which is the general term of a converging series. 

Passing to the limit as n — > +00, r — > so that 2rn — > cr, one finds 
(by dominated convergence in the series) that 

(f r M (t,x,v,{c})) ^ e-^f in (x -tv,v) 

nt p2tt 

+<je-° t / / f m (x - Tl v -(t- r^R^v, R\fc]v) sin f *&dn 
Jo Jo 



+Y t ° je ~ ot f f 

J0<T,<...<TT<t J\C 



X 



-Jj\r k R 



k=l 



i>2 

k-l 

j=i 



0<-r j <...<Ti<t J[0,2vr]i 



sin §- . . . sin § 



5> 

j=i 



^j-drx . . . d,Tj 



4 ••• 4 "'i ,,,u "3i 

which is the Duhamel series giving the solution of the Lorentz kinetic 
equation. 

Hence, we have proved that 

{f™{t,x,v, •)) — > f(t,x,v) uniformly on bounded sets as r — > + , 

where / is the solution of the Lorentz kinetic equation. One can check 
by a straightforward computation that the Lorentz collision integral 
satisfies the property 



C{4>){v)dv = for each G ^(S 1 ) . 



Integrating both sides of the Lorentz kinetic equation in the variables 
(t, x, v) over [0, t] x R 2 x S 1 shows that the solution / of that equation 
satisfies 



f(t, x, v)dxdv 



R2XS 1 



f m (x,v)dxdv 



R2XS 1 



for each t > 0. 

On the other hand, the billiard flow (X, V){t,-,-, {c}) obviously leaves 
the uniform measure dxdv on R 2 x S 1 (i.e. the particle number) in- 
variant, so that, for each t > and each r > 0, 



SL 



R 2 xS x 



f r (t, x, v, {c})dxdv 



R 2 xS x 



f m {x,v)dxdv . 



We therefore deduce from Fatou's lemma that 

(f?) — > in L X (R 2 x S 1 ) uniformly on bounded t-sets 

(/ r M ) — > / in L X (R 2 x S 1 ) uniformly on bounded t-sets 

which concludes our sketch of the proof of Gallavotti's theorem. 



□ 
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For a complete proof, we refer the interested reader to [HJ [20] . 

Some remarks are in order before leaving Gallavotti's setting for the 
Lorentz gas with the Poisson distribution of obstacles. 

Assuming no external force field as done everywhere in the present 
paper is not as inocuous as it may seem. For instance, in the case of 
Poisson distributed holes — i.e. purely absorbing obstacles, so that 
particles falling into the holes disappear from the system forever - 
the presence of an external force may introduce memory effects in the 
Boltzmann-Grad limit, as observed by L. Desvillettes and V. Ricci |15j . 

Another remark is about the method of proof itself. One has ob- 
tained the Lorentz kinetic equation after having obtained an explicit 
formula for the solution of that equation. In other words, the equation 
is deduced from the solution — which is a somewhat unusual situation 
in mathematics. However, the same is true of Lanford's derivation of 
the Boltzmann equation [2H], as well as of the derivation of several 
other models in nonequilibrium statistical mechanics. For an interest- 
ing comment on this issue, see [13], on p. 75. 

3. Santalo's formula 
for the geometric mean free path 

From now on, we shall abandon the random case and concentrate 
our efforts on the periodic Lorentz gas. 

Our first task is to define the Boltzmann-Grad scaling for periodic 
systems of spherical obstacles. In the Poisson case defined above, things 
were relatively easy: in space dimension 2, the Boltzmann-Grad scaling 
was defined by the prescription that the number of obstacles per unit 
volume tends to infinity while the obstacle radius tends to in such a 
way that 

# obstacles per unit volume x obstacle radius — > a > . 

The product above has an interesting geometric meaning even with- 
out assuming a Poisson distribution for the obstacle centers, which we 
shall briefly discuss before going further in our analysis of the periodic 
Lorentz gas. 

Perhaps the most important scaling parameter in all kinetic models 
is the mean free path. This is by no means a trivial notion, as will be 
seen below. As suggested by the name itself, any notion of mean free 
path must involve first the notion of free path length, and then some 
appropriate probability measure under which the free path length is 
averaged. 
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Figure 6. The periodic billiard table 

For simplicity, the only periodic distribution of obstacles considered 
below is the set of balls of radius r centered at the vertices of a unit 
cubic lattice in the D- dimensional Euclidian space. 

Correspondingly, for each r G (0 ,|), we define the domain left free 
for particle motion, also called the "billiard table" as 

Z T = {x G R D | dist(x, Z D ) > r} . 

Defining the free path length in the billiard table Z r is easy: the free 
path length starting from x G Z r in the direction v G S D_1 is 

r r (x, v) = min{t > | x + tv G dZ r } . 

Obviously, for each v G S D_1 the free path length r r (-,v) in the 
direction v can be extended continuously to 

{x G dZ r | v ■ n x 7^ 0} , 

where n x denotes the unit normal vector to dZ r at the point x G dZ r 
pointing towards Z r . 

With this definition, the mean free path is the quantity defined as 

Mean Free Path = (r r ) , 

where the notation (•) designates the average under some appropriate 
probability measure on Z r x S D_1 . 

A first ambiguity in the notion of mean free path comes from the fact 
that there are two fairly natural probability measures for the Lorentz 
gas. 
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\(x,v) 



Figure 7. The free path length 

The first one is the uniform probability measure on Z r /Z D x S D_1 

dxdv 

dfj, r (x, v) 



\z r /z D \ (s^- 1 ! 

that is invariant under the billiard flow — the notation |S D_1 | desig- 
nates the D — 1-dimensional uniform measure of the unit sphere S ^ 1 . 
This measure is obviously invariant under the billiard flow 

(X r , V r ^j(t, •, •) ; Z r x S^* ^ — ► Z r x 

defined by 

~ ^ r whenever X(t) dZ r 
while 

X r (t + ) = X r (t~) =: X r (t) if X(t ± ) e dZ r , 

V r (t+) =n[n Xr(t) ]v r {t-) 

with TZ[n]v = v — 2v-nn denoting the reflection with respect to the 
hyperplane (Rn)- 1 . 

The second such probability measure is the invariant measure of the 
billiard map 

, , . v-n x dS(x)dv 

dh>Jx,v) = : — : 7— — 

v ■ n x dxdv-meas (T\ /Z u ) 

where n x is the unit inward normal at x G dZ r , while dS(x) is the 
D — 1-dimensional surface element on dZ r , and 

r; := {(x,v) e dZ r x S°- l \vn x > 0}. 

The billiard map B r is the map 

3 (x, v) i-> B r {x, v) := {X r , V r )(T r (x, v); x, v) E T r + , 

which obviously passes to the quotient modulo Z D -translations: 

B r : r + /z D -> r;/z D . 

In other words, given the position x and the velocity v of a parti- 
cle immediatly after its first collision with an obstacle, the sequence 
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(Br(x,v))n>Q is the sequence of all collision points and post-collision 
velocities on that particle's trajectory. 

With the material above, we can define a first, very natural notion 
of mean free path, by setting 

1 N ^ 

Mean Free Path = lim — r r (B^(x,v)) . 

Notice that, for u r -a.e. (x,v) G r+/Z D , the right hand side of the 
equality above is well-defined by the Birkhoff ergodic theorem. If the 
billiard map B r is ergodic for the measure z/ r , one has 

1 N ^ f 

lim — ) T r (Br(x,v)) = / T r du r , 

^+°° N Jr\m D 



N 

k=0 



for u r -a.e. (x,v) G T^/Z . 

Now, a very general formula for computing the right-hand side of 
the above equality was found by the great Spanish mathematician L. 
A. Santalo in 1942. In fact, Santalo's argument applies to situations 
that are considerably more general, involving for instance curved trajec- 
tories instead of straight line segments, or obstacle distributions other 
than periodic. The reader interested in these questions is referred to 
Santalo's original article 

Here is 



1 - 


B 


D 




\B°- 


-i 


r D-l 



Santalo's formula for the geometric mean free path 
One finds that 

£ r := / T r (x,v)du r (x,v) 
Jr r + /z D 

where B D is the unit ball of R/° and \B D \ its D-dimensional Lebesgue 
measure. 

In fact, one has the following slightly more general 

Lemma 3.1 (H.S. Dumas, L. Dumas, F. Golse [T7]). For f e C 1 (R+) 
such that /(0) = 0, one has 



f(r r (x,v))v ■ n x dS{x)dv = // f{r r {x,v))dxdv . 

r r + /z D J J(z r /z D )xs D - 1 

Santalo's formula is obtained by setting f(z) = z in the identity 
above, and expressing both integrals in terms of the normalized mea- 
sures v r and fx r . 
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Figure 8. Luis Antonio Santalo Sors (1911-2001) 



Proof. For each (x,v) G Z r x S D_1 one has 

T r (x + tv, v) = T r (x, V ) — t , 

so that 

d 

— T r {X + tv , V J = — 1 . 

at 

Hence r r (x,v) solves the transport equation 

v ■ V x r r (x, v ) = — 1 , x G Z r , v G S D_1 , 
T r (x, v ) = , a; G <9Z r , t> • n x < . 

Since / G C 1 (R + ) and /(0) = 0, one has 

f « • V,/(r r (x, «)) = -f(r r (x, v)), x G Z r , v G S^ 1 , 
\ f(r r (x, v)) = , x G <9Z r , u • n x < 
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Integrating both sides of the equality above, and applying Green's for- 
mula shows that 



f'(r r (x, v))dxdv 

(Z r /Z D )xS D -i 



v ■ V x (f(r r (x,v)))dxdv 

(Z r /Z D )XS D - 1 

f(r r (x, v))v ■ n x dS(x)dv 

{dZ r /Z D )xS D - 1 

- beware the unusual sign in the right-hand side of the second equality 
above, coming from the orientation of the unit normal n x , which is 
pointing towards Z r . □ 

With the help of Santalo's formula, we define the Boltzmann-Grad 
limit for the Lorentz gas with periodic as well as random distribution 
of obstacles as follows: 



Boltzmann-Grad scaling 



The Boltzmann-Grad scaling for the periodic Lorentz gas in space 
dimension D corresponds with the following choice of parameters: 

distance between neighboring lattice points = e , 

obstacle radius = r < 1, 

mean free path = £ r — > — > . 

a 

Santalo's formula indicates that one should have 



d i a 



D-l 



r ~ ce - 1 with c = ^ jgD-ij J as £ — > . 

Therefore, given an initial particle density f m G C C (R D x S D_1 ), we 
define f r to be 

WW) = r (r D -'X r ,V r (-^; ^,,)) 

where (X r , V r ) is the billiard flow in Z r with specular reflection on dZ r . 

Notice that this formula defines f r for x G Z r only, as the particle 
density should remain for all time in the spatial domain occupied by 
the obstacles. As explained in the previous section, this is a set whose 
measure vanishes in the Boltzmann-Grad limit, and we shall always 
implicitly extend the function f r defined above by for x Z r . 
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Since f m is a bounded function on Z r x S D_1 , the family f r de- 
fined above is a bounded family of /^(R 15 x S* " 1 ). By the Banach- 
Alaoglu theorem, this family is therefore relatively compact for the 
weak-* topology of L°°(R + x R D x S^ 1 ). 

Problem: to find an equation governing the L°° weak-* limit points 
of the scaled number density f r as r — > + . 

In the sequel, we shall describe the answer to this question in the 
2-dimensional case (D = 2.) 

4. Estimates for the distribution of free-path lengths 

In the proof of Gallavotti's theorem for the case of a Poisson dis- 
tribution of obstacles in space dimension D = 2, the probability that 
a strip of width 2r and length t does not meet any obstacle is e~ 2nrt , 
where n is the parameter of the Poisson distribution — i.e. the average 
number of obstacles per unit surface. 

This accounts for the loss term 

f in {x-tv,v)e- (Tt 

in the Duhamel series for the solution of the Lorentz kinetic equation, 
or of the term —erf on the right-hand side of that equation written in 
the form 

(d t + v ■ V x )f = -crf + a f(t, x, R(J3)v) sin ff . 

JO 

Things are fundamentally different in the periodic case. To begin 
with, there are infinite strips included in the billiard table Z r which 
never meet any obstacle. 

The contribution of the 1-particle density leading to the loss term 
in the Lorentz kinetic equation is, in the notation of the proof of 
Gallavotti's theorem 

f m (x -tv,v)l t<Tl{xMc}) . 
The analogous term in the periodic case is 

f m (x - tV, v)l t<r D-i Tr ( x / r _ v <) 

where r r (x,v) is the free-path length in the periodic billiard table Z r 
starting from x G Z r in the direction v G S 1 . 

Passing to the L°° weak-* limit as r — > reduces to finding 

liml t<r r>-.i Tr(a ./ r _„) in w* - L°°(R 2 x S 1 ) 

- possibly after extracting a subsequence r n j 0. As we shall see 
below, this involves the distribution of r r under the probability measure 
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Figure 9. Open strips in the periodic billiard table that 
never meet any obstacle 

\x r introduced in the discussion of Santalo's formula — i.e. assuming 
the initial position x and direction v to be independent and uniformly 
distributed on (R D /Z D ) x S^ 1 . 

We define the (scaled) distribution under /i r of free path lengths r r 
to be 

$ r (t) := (jtr({(x,v) e {Z r /Z D ) x S^lrrC^v) >t/r D - 1 }). 

Notice the scaling t i— > t/r D ~ l in this definition. In space dimension 
D, Santalo's formula shows that 

// T r (x,v)dv r (x,v) ~ |Bi^T| rl_D > 

JJr+/zo 1 1 

and this suggests that the free path length r r is a quantity of the order 
of 1/r . (In fact, this argument is not entirely convincing, as we 
shall see below.) 

In any case, with this definition of the distribution of free path 
lengths under fi r , one arrives at the following estimate. 

Theorem 4.1 (Bourgain-Golse-Wennberg [6JI25]). In space dimension 
D > 2, there exists < Co < C' D such that 

< <& r (t) < — — whenever t > 1 and < r < \ . 

t ~ y 1 ~ t 2 

The lower bound and the upper bound in this theorem are obtained 
by very different means. 
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The upper bound follows from a Fourier series argument which is 
reminiscent of Siegel's prood of the classical Minkowski convex body 
theorem (see [39| 136].) 

The lower bound, on the other hand, is obtained by working in phys- 
ical space. Specifically, one uses a channel technique, introduced inde- 
pendently by P. Bleher [2] for the diffusive scaling. 

This lower bound alone has an important consequence: 

Corollary 4.2. For each r > 0, the average of the free path length 
(mean free path) under the probability measure fi r is infinite: 

/ T r (x, v)dfj, r (x, v ) = +00 . 

J(Z r /Z D )xS D - 1 

Proof. Indeed, since $ r is the distribution of r r under /i r , one has 



L 



oo c 



r r (x, v)d/i r (x, v) = I $> r (t)dt > I — —dt = +oo . 



{Z r /Z D )xS D - 1 ' JO Jl t 

□ 

Recall that the average of the free path length unded the "other" 
natural probability measure v r is precisely Santalo's formula for the 
mean free path: 



1 - 


B 


D\ 






-l 


r D-l 



r r (x, v)dv r (x, v) 

r+/z D 

One might wonder why averaging the free path length r r under the 
measures u r and fi r actually gives two so different results. 

First observe that Santalo's formula gives the mean free path under 
the probability measure v r concentrated on the surface of the obstacles, 
and is therefore irrelevant for particles that have not yet encountered 
an obstacle. 

Besides, by using the lemma that implies Santalo's formula with 

f(z) = \z 2 , one has 

r r (x, v)d/JL r (x, v) = — I \t t (x, v ) 2 du r (x, V ) . 

{Z r /Z D )xS D - 1 Jr+/z D 

Whenever the components vi,...,Vd are independent over Q, the 
linear flow in the direction v is topologically transitive and ergodic on 
the .D-torus, so that r r (x, v) < +oo for each r > and x G R D . On the 
other hand, r r (x, v) = +oo for some x E Z r (the periodic billiard table) 
whenever v belongs to some specific class of unit vectors whose compo- 
nents are rationally dependent, a class that becomes dense in S D_1 as 
r — ► + . Thus, r r is strongly oscillating (finite for irrational directions, 
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possibly infinite for a class of rational directions that becomes dense 
as r — > + ), and this explains why r r doesn't have a second moment 
under u r . 

Proof of the Bourgain-Golse-Wennberg lower bound 

We shall restrict our attention to the case of space dimension D = 2. 

As mentionned above, there are infinite open strips included in Z r 
- i.e. never meeting any obstacle. Call a channel any such open strip 
of maximum width, and let C r be the set of all channels included in Z r . 

If S G C r and x G S, define Ts(x,v) the exit time from the channel 
starting from x in the direction v, defined as 

t s (x, v) = inf {£ > | x + tv G OS} , (x, v) G S x S 1 . 

Obviously, any particle starting from x in the channel S in the direction 
v must exit £ before it hits an obstacle (since no obstacle intersects S). 
Therefore 

T r (x, v) > sup{rs(x,v) | S G C r s.t. 

so that 

$ r (t) >ii r ( \J{(x,v) G (5/Z 2 ) x S 1 |r s (a;,7;) > t/r} 
\sec r 

This observation suggests that one should carefully study the set of 
channels C r . 

Step 1: description of C r . Given u G S 1 , we define 

C r (u) := {channels of direction u in C r } ; 

We begin with a lemma which describes the structure of C r {uj). 

Lemma 4.3. Let r G [0, \) and uo G S 1 . T/ien 
if S £ C r (oo), then 

C r (u) := {S + k\ k G Z 2 }; 

^ ifC r (uj) ^ 0, then 

(p, q) 

uj = 

VP 2 + q 2 

with 

(v, q) G Z 2 \ {(0, 0)} such that g.c.d.(p, q) — 1 and \/ p 2 + q 2 < — . 
We henceforth denote by A r the set of all such u G S 1 . T/ien 




PERIODIC LORENTZ GAS 



25 



2r J 
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L 1 

FIGURE 10. A channel of direction uj = ^(2,1); min- 
imal distance d between lines L and V of direction uj 
through lattice points 



3) for uj e A r , the elements of C r (u) are open strips of width 

1 



w(uj, r) 



yV + q 2 



2r. 



Proof of the Lemma. Statement 1) is obvious. 

As for statement 2), if L is an infinite line of direction w G S 1 such 
that U2/UJ1 is irrational, then L/Z 2 is an orbit of a linear flow on T 2 
with irrational slope u^/^i- Therefore L/Z 2 is dense in T 2 so that L 
cannot be included in Z r . 

Assume that 

(p,q) 



with (p,q) 6Z \ {(0,0)} coprime, 



and let L, V be two infinite lines with direction uj, with equations 
qx — py = a and qx — py = a' respectively. 

Obviously 

dist(L,L') 



\a — a 



\Zp 2 + ( i 2 

If L U L' is the boundary of a channel of direction 



26 F. GOLSE 

included in R 2 \Z 2 - - i.e. of an element of Cq(u), then L and L' 
intersect Z 2 so that 

a, a' G pZ + qZ = Z 

- the equality above following from the assumption that p and q are 
coprime. 

Since dist(L, L') > is minimal, then \a — a'\ = 1, so that 
dist(L,L') ! 



vV + <? 2 

Likewise, if L U L' = with S G C r , then L and L' are parallel 
infinite lines tangent to dZ r , and the minimal distance between any 
such distinct lines is 

dirt(L,L') = -— L=-2r. 
VP + Q 

This entails 2) and 3). □ 



Step 2: the exit time from a channel. Let ui = ) p ' <? - > e »4 r and let 
VP 2 +q 2 

S G C r (c<j). Cut 5 into three parallel strips of equal width and call S 

the middle one. For each t > 1 define 



= 9(uj,r,t) := arcsin 



3t 



Lemma 4.4. //x G 5 and v G 0]c<j, i?[0]c<j) ; where R[6] designates 
the rotation of an angle 6, then 

ts{x,v) >t/r. 

Moreover 

l2 r ((S/Z 2 ) x (R[-6]uj,R[e]uj)) = lw(u,r)6(u,r,t). 

The proof of this lemma is perhaps best explained by considering 
Figure 11. 

Step 3: putting all channels together. Recall that we need to estimate 

Hr ( |J {(x,v) G (S/Z 2 ) x S 1 I t s (x,v) > t/r}) . 
\sec r J 



Pick 



A r 3 U = —7== ± -^=0 = £ A r 
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Figure 11. Exit time from the middle third S of an 
infinite strip S of width w 



Observe that 



sin w, uS 



\pq'-p'q\ 



> 



y/p 2 +q 2 y/p' 2 +q' 2 ~ ^P 2 +q 2 y/p' 2 +q' 2 



> max ( 



2r 



> sin 8(u, r, t) + sin 8(u', r, t) 



_ y/p 2 +q 2 ' y/p' 2 +q' 2 

> sm(9(u,r,t) + 9(u/,r,t)) 

whenever t > 1. 

Then, whenever S G C r (u) and S' G C r (u') 

(S x (R[-6]u},R[0]u))) n (5' x (i^V^V))) = 

with # = 9(u,r,t), 9' = 9'(u',r,t) and R[9] =the rotation of an angle 
9. 

Moreover, if uj = ^ p ' q ^ G A r then 
yp 2 +q 2 



while 



\S/Z 2 \ = \w(cu,r)^p 2 + q 2 . 



#{S/Z 2 \S eC r {uo)} = 1. 



Conclusion: Therefore, whenever t > 1 



(J (S/Z 2 ) x (R[-0]u,R[9]u) 

C |J {(x, v) G (5/Z 2 ) x S 1 1 r s {x, v) > t/r} . 



sec 



sec,- 
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A channel modulo Z 2 

and the left-hand side is a disjoint union. Hence 

fir ( |J {(x, v) e (S/Z 2 ) x S 1 1 t s (x, v) > t/r} 
\s<=c r 

> J2 Vr((S/Z 2 ) x (R[-9}u,R[6}u)) 
= lw(u,r)^/p 2 + q 2 -29(u,r,t) 

rw(u, r) 
3t 

rw(u>, r) 2 



g.c.d.(p,q)=X 
p 2 +q 2 <l/4r 2 



/J |a/p 2 + q 2 w(u, r) arcsin I 

d.(p,q) = l \ 



g.c.d.(p,q) = 
p 2 +q 2 <l/4r 2 



> E Iv^r 



g.c.d.(p,ij) = l 
p 2 +<j 2 <l/4r 2 



3t 



Now 



\J p 2 + q 2 < l/4r =>• uj(a>, r) 



1 - 2r > 1 
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FIGURE 12. Black lines issued from the origin terminate 
at integer points with coprime coordinates; red lines ter- 
minate at integer points whose coordinates are not co- 
prime 



so that, eventually 



g.c.d.(p,q) = l 
p 2 +<2 2 <l/16r 2 



: ra(w, r)' 
3t 



> 



18* 



E 



g.c.d.(p,q) = l 
p 2 + q 2 <l/16r 2 



Vf 2 + 



This gives the desired conclusion since 

1 



E 



g.c.d.(p,q) = l 
p 2 +q 2 <l/ldr 2 



Ar^Jp 2 + q 2 



E 



7T 



p 2 +q 2 <l/16r 2 



16r 



2 ' 



The first equality is proved as follows: the term 

1 



4r a/p 2 + q 2 



is the number of integer points on the segment of length l/4r in the 
direction (p,q) with (p, q) G Z 2 such that g.c.d.(p,q) = 1. 
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The Bourgain-Golse-Wennberg theorem raises the question, of whe- 
ther $ r (t) ~ C/t in some sense as r — > + and t — > +00. Given the 
very different nature of the arguments used to establish the upper and 
the lower bounds in that theorem, this is a highly nontrivial problem, 
whose answer seems to be known only in space dimension D = 2 so far. 
We shall return to this question later, and see that the 2-dimensional 
situation is amenable to a class of very specific techniques based on 
continued fractions, that can be used to encode particle trajectories of 
the periodic Lorentz gas. 

A first answer to this question, in space dimension D = 2, is given 
by the following 

Theorem 4.5 (Caglioti-Golse [9]). Assume D = 2 and define, for each 
r C S . 

(j) r (t\v) = fx r ({x E Z r /Z 2 I r r (x, v) > t/r} , t > . 
Then there exists $ : R + — -> R + such that 

1 r 1/4 , Ar 



ft dr 

I (f>r(t,v) > $(t) a.e. inveS 1 

Je r 



I lne 

in the limit as e — > + . Moreover 

$(t) ~ ast^ +00 . 

Shortly after [5] appeared, F. Boca and A. Zaharescu improved our 
method and managed to compute explicitly for each t > 0. One 
should keep in mind that their formula had been conjectured earlier by 
P. Dahlqvist [14J , on the basis of a formal computation. 

Theorem 4.6 (Boca-Zaharescu 0). For each t > 

/CO 
(s-t)g(s)ds 

in the limit as r — > + , where 



_ 1 s G [0, 1] 

^ (S) "^ Xi i + 2(1- i) 2 lrx(l - i) - 1 11 - f | 2 lxx|l - f I *e(l,oo) 
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In the sequel, we shall return to the continued and Farey fractions 
techniques used in the proofs of these two results, and generalize them. 
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Figure 13. Graph of $(t) (right) and g(t) = $"(t) (left) 

5. A NEGATIVE RESULT FOR THE BOLTZMANN-GRAD LIMIT 
OF THE PERIODIC LORENTZ GAS 

The material at our disposal so far provides us with a first answer - 
albeit a negative one — to the problem of determining the Boltzmann- 
Grad limit of the periodic Lorentz gas. 

For simplicity, we consider the case of a Lorentz gas enclosed in a 
periodic box T D = R D /Z D of unit side. The distance between neigh- 
boring obstacles is supposed to be e ^ 1 with < e — 1/n, for n G N 
and n > 2 so that e < 1/2, while the obstacle radius is e D < \s D ~ 1 - 
so that obstacles never overlap. Define 

Y e = {xe T D |dist(ar,e iJ - 1 Z IJ ) > e D } = s D -\Z £ /Z D ) . 

For each f in G C(T D x S D ^), let f £ be the solution of 

dtfe + v- V x f £ = 7 (x,v)eY £ xS D - 1 

f £ (t, x, v) = f e (t, x,K[n x ]v) , (x, v) G dY £ x S ' 1 

r I ?in 

Je\t=0 ~ J ' 

where n x is unit normal vector to dY £ at the point x, pointing towards 
the interior of Y £ . 

By the method of characteristics 

f e (t,x,v) = f m {e D - 1 X £ (-^; ^v) ; V £ (-^; ^=r,v)) 

where (X £ , V £ ) is the billiard flow in Z £ . 

The main result in this section is the following 

Theorem 5.1 (Golse (2111211). There exist imtial data f n = f in ( x ) e 
C(T D ) such that no subsequence of f £ converges for the weak-* topology 
of L°°(R + xT D xS ^ 1 ) to the solution f of a linear Boltzmann equation 



32 



F. GOLSE 



of the form 



(dt + v ■ V x )f(t, x,v) = a \ p(v,v')(f(t,x,v') - f(t,x,v))dv' 



Js - 1 




! p(v,v')dv' = / p(v',v)dv' = l a.e. in v 6 S 

SO- 1 Jsd- 1 



This theorem has the following important — and perhaps surpris- 
ing — consequence: the Lorentz kinetic equation cannot govern the 
Boltzmann-Grad limit of the particle density in the case of a periodic 
distribution of obstacles. 

Proof. The proof of the negative result above involves two different 
arguments: 

a) the existence of a spectral gap for any linear Boltzmann equation, 
and 

b) the lower bound for the distribution of free path lengths in the 
Bourgain-Golse-Wennberg theorem. 

Step 1: Spectral gap for the linear Boltzmann equation: 



With cr > and p as above, consider the unbounded operator A on 
L 2 (T D xS w ) defined by 



D(A) = {0 G L 2 (T D x S D - X ) | v ■ V x (j> G L 2 (T D x S^ 1 )} . 

Then 



Theorem 5.2 (Ukai-Point-Ghidouche |45j). There exists positive con- 
stants C and 7 such that 



||e* (p~ (0)||l2(tb x sd-i) < Ce ^H^Il^xS - 1 ) , *>0, 
for each <\> G L 2 (T D x S ' 1 ), where 



J JT D xS D ~ 1 

Taking this theorem for granted, we proceed to the next step in the 
proof, leading to an explicit lower bound for the particle density. 

Step 2: Comparison with the case of absorbing obstacles 




with domain 




4>(x, v)dxdv . 
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Assume that f in = f in (x) > on T D . Then 

f e (t,X,v) > g £ (t,X,v) = f m (x - tv)l Ye (x)l e D-i Te ( x/E D-i jV)>t . 

Indeed, g is the density of particles with the same initial data as /, but 
assuming that each particle disappear when colliding with an obstacle 
instead of being reflected. 
Then 

l Ye (x) -> 1 a.e. on T D and \l Ye ( x )\ < 1 
while, after extracting a subsequence if needed, 

h^-W e (x/ e ^,v)>t ~* *(t,v) in L°°(R + x T D x S^ 1 ) weak-* . 

Therefore, if / is a weak-* limit point of f £ in L°°(R + x T D x S ' 1 ) 
as e — > 

/(*, ^) > /* "(z - tv)^(t, v) for a.e. (t, x, u) . 

Step 3: using the lower bound on the distribution of r r 
Denoting by dv the uniform probability measure on S D_1 

f(t, x, v) 2 dxdv 



\s°-i\ 



T D xS D-l 



> js^u [ [ f n {x - tv) 2 *(t, vfdxdv 

= / f n {y) 2 dy^hu [ *(t,v) 2 dv 



By the Bourgain-Golse-Wennberg lower bound on the distribution $ 
of free path lengths 

On the other hand, by the spectral gap estimate, if / is a solution of 
the linear Boltzmann equation, one has 

\\f(t,;-)\\ LH T°x S °-i)< [ f in (y)dy + Ce-*\\f in \\ LHTD) 

Jt d 



so that 



for each t > 1. 



Cd_ < ll/ m |lLi(T^) Ce _ 7t 
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Step 4: choice of initial data 

Pick p to be a bump function supported near x = and such that 

p{x) 2 dx = 1 . 

Take f m to be x i— > A D / 2 p(Ax) periodicized, so that 

I f n {xfdx = 1, while / f in (y)dy = A~ D/2 I p{x)dx . 
Jt d Jt d J 

For such initial data, the inequality above becomes 

Cd 
t 

Conclude by choosing A so that 



< \- D/2 J p{x)dx + Ce~ 7t . 



\- D/2 J p(x)dx < sup (^j- - Ce^J > . 



□ 



Remarks: 

1) The same result (with the same proof) holds for any smooth obstacle 
shape included in a shell 

{x G R D | Ce D < dist(x,£ Z? - 1 Z D ) < C'e D } . 

2) The same result (with same proof) holds if the specular reflection 
boundary condition is replaced by more general boundary conditions, 
such as absorption (partial or complete) of the particles at the bound- 
ary of the obstacles, diffuse reflection, or any convex combination of 
specular and diffuse reflection — in the classical kinetic theory of gases, 
such boundary conditons are known as "accomodation boundary con- 
ditions" . 

3) But introducing even the smallest amount of stochasticity in any 
periodic configuration of obstacles can again lead to a Boltzmann-Grad 
limit that is described by the Lorentz kinetic model. 

Example. (Wennberg-Ricci [37]) In space dimension 2, take obsta- 
cles that are disks of radius r centered at the vertices of the lattice 
r 1 /( 2 -';)Z 2 ) assuming that < r\ < 1. SantakVs formula suggests that 
the free-path lengths scale like r^ 2- *^ — > 0. 

Suppose the obstacles are removed independently with large prob- 
ability — specifically, with probability p = \ — r r, ^ 2 ~ v \ In that case, the 
Lorentz kinetic equation governs the 1-particle density in the Boltzmann- 
Grad limit as r — ► + . 

Having explained why neither the Lorentz kinetic equation nor any 
linear Boltzmann equation can govern the Boltzmann-Grad limit of the 
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periodic Lorentz gas, in the remaining part of these notes, we build the 
necessary material used in the description of that limit. 

6. Coding particle trajectories with continued fractions 

With the Bourgain-Golse-Wennberg lower bound for the distribution 
of free path lengths in the periodic Lorentz gas, we have seen that the 
1-particle phase space density is bounded below by a quantity that is 
incompatible with the spectral gap of any linear Boltzmann equation 
— in particular with the Lorentz kinetic equation. 

In order to further analyze the Boltzmann-Grad limit of the periodic 
Lorentz gas, we cannot content ourselves with even more refined esti- 
mates on the distribution of free path lengths, but we need a convenient 
way to encode particle trajectories. 

More precisely, the two following problems must be answered some- 
how: 

First problem: for a particle leaving the surface of an obstacle in a given 
direction, to find the position of its next collision with an obstacle; 

Second problem: average — in some sense to be defined — in order to 
eliminate the direction dependence. 

From now on, our discussion is limited to the case of spatial dimen- 
sion D = 2, as we shall use continued fractions, a tool particularly well 
adapted to understanding the rational approximation of real numbers. 
Treating the case of a space dimension D > 2 along the same lines 
would require a better understanding of simultaneous rational approx- 
imation of D — 1 real numbers (by D — 1 rational numbers with the 
same denominator), a notoriously more difficult problem. 

We first introduce some basic geometrical objects used in coding 
particle trajectories. 

The first such object is the notion of impact parameter. 

For a particle with velocity v G S 1 located at the position x on the 
surface of an obstacle (disk of radius r), we define its impact parameter 
h r (x,v) by the formula 

h r (x,v) = sin(n^py) . 

In other words, the absolute value of the impact parameter h r (x,v) is 
the distance of the center of the obstacle to the infinite line of direction 
v passing through x . 
Obviously 

h r (x, 1Z[n x ]v) = h r (x,v) 
where we recall the notation TZ[n]v = v — 2v ■ nn. 
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Figure 14. The impact parameter h corresponding 
with the collision point x at the surface of an obstacle, 
and a direction v 




FIGURE 15. The transfer map 



The next important object in computing particle trajectories in the 
Lorentz gas is the transfer map. 

For a particle leaving the surface of an obstacle in the direction v 
and with impact parameter h', define 

ji (ft v ) — ( s m w jth / S = rX distance ^° ^ ne nex ^ collision point 
r ^ ' ' ^ ' ' 1 h = impact parameter at the next collision 

Particle trajectories in the Lorentz gas are completely determined by 
the transfer map T r and its iterates. 

Therefore, a first step in finding the Boltzmann-Grad limit of the 
periodic, 2-dimensional Lorentz gas, is to compute the limit of T r as 
r — ► + , in some sense that will be explained later. 
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FIGURE 16. Three types of orbits: the blue orbit is the 
shortest, the red one is the longest, while the green one 
is of the intermediate length. The black segment removed 
is orthogonal to the direction of the trajectories. 

At first sight, this seems to be a desperately hard problem to solve, 
as particle trajectories in the periodic Lorentz gas depend on their 
directions and the obstacle radius in the strongest possible way. For- 
tunately, there is an interesting property of rational approximation on 
the real line that greatly reduces the complexity of this problem. 

The 3-length theorem 

Question (R. Thom, 1989): on a flat 2-torus with a disk removed, 
consider a linear flow with irrational slope. What is the longest orbit? 

Theorem 6.1 (Blank-Krikorian [lj). On a flat 2-torus with a segment 
removed, consider a linear flow with irrational slope < a < 1. The 
orbits of this flow have at most 3 different lengths — exceptionally 2, 
but generically 3. Moreover, in the generic case where these orbits have 
exactly 3 different lengths, the length of the longest orbit is the sum of 
the two other lengths. 

These lengths are expressed in terms of the continued fraction ex- 
pansion of the slope a. 
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Figure 17. The 3-term partition. The shortest orbits 
are collected in the blue strip, the longest orbits in the red 
strip, while the orbits of intermediate length are collected 
in the green strip. 



Together with E. Caglioti in [9], we proposed the idea of using the 
Blank-Krikorian 3-length theorem to analyze particle paths in the 2- 
dimensional periodic Lorentz gas. 

More precisely, orbits with the same lengths in the Blank-Krikorian 
theorem define a 3-term partition of the flat 2-torus into parallel strips, 
whose lengths and widths are computed exactly in terms of the con- 
tinued fraction expansion of the slope (see Figure 17q) 

The collision pattern for particles leaving the surface of one obstacle 
- and therefore the transfer map — can be explicitly determined in 
this way, for a.e. direction v G S 1 . 

In fact, there is a classical result known as the 3-length theorem, 
which is related to Blank-Krikorian's. Whereas the Blank-Krikorian 
theorem considers a linear flow with irrational slope on the flat 2-torus, 

Figures 16 and 17 are taken from a conference by E. Caglioti at the Centre 
International de Rencontres Mathematiques, Marseilles, February 18-22, 2008. 




the classical 3-length theorem is a statement about rotations of an 
irrational angle — i.e. about sections of the linear flow with irrational 
slope. 

Theorem 6.2 (3-length theorem). Let a <E (0, 1) \ Q and N > 1. The 

sequence 

{na | < n < N} 

defines iV + 1 intervals on the circle of unit length ~ R/Z. The lengths 
of these intervals take at most 3 different values. 

This striking result was conjectured by H. Steinhaus, and proved in 
1957 independently by P. Erdos, G. Hajos, J. Suranyi, N. Swieczkowski, 
P. Sziisz — reported in [12], and by Vera Sos [JT]. 

As we shall see, the 3-length theorem (in either form) is the key to 
encoding particle paths in the 2-dimensional Lorentz gas. We shall need 
explicitly the formulas giving the lengths and widths of the 3 strips in 
the partition of the flat 2-torus defined by the Blank-Krikorian theorem. 
As this is based on the continued fraction expansion of the slope of the 
linear flow considered in the Blank-Krikorian theorem, we first recall 
some basic facts about continued fractions. An excellent reference for 
more information on this subject is [25] . 

Continued fractions 
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Assume < i>2 < v\ and set a = i^/fi, and consider the continued 
fraction expansion of a: 

r i 1 

a = [0; a , ai, a 2 , . . .J = j . 

a + 



ai + ... 

Define the sequences of convergents (p n , q n )n>o — meaning that 

Pn+2 r n 1 \ o 

= [0; ao, ■ ■ ■ , a n \ , n>2 

Qn+2 

— by the recursion formulas 

p n+ i = a n p n + p n _i , po = 1 , pi = , 
q n+ i = a n q n + q n _ x q = , ft = 1 . 

Finally, let d„ denote the sequence of errors 

d n = \q n OL ~ p n \ = {-l) n ~ l {q n a -Pn), U > , 

so that 

d n +i — —dndn + d n _\ , do — 1 , d\ — a . 
The sequence d n is decreasing and converges to 0, at least exponen- 
tially fast. (In fact, the irrational number for which the rational approx- 
imation by continued fractions is the slowest is the one for which the 
sequence of denominators q n have the slowest growth, i.e. the golden 
mean 

l VE-i 



61 = [0; 1,1,...] 



1 + 



1 + 

The sequence of errors associated with 6 satisfies d n+ i = —d n + 
for each n > 1 with d — 1 and d 1 = 9, so that d n = 9 n for each n > 0.) 
By induction, one verifies that 

q n d n+1 + q n+1 d n = 1 , n > . 

Notation: we write p n (a) , q n (a) , d n (a) to indicate the dependence of 
these quantities in a. 
Collision patterns 

The Blank- Krikorian 3-length theorem has the following consequence, 
of fundamental importance in our analysis. 

Any particle leaving the surface of one obstacle in some irrational 
direction v will next collide with one of at most 3 — exceptionally 2 - 
other obstacles. 

Any such collision pattern involving the 3 obstacles seen by the de- 
parting particle in the direction of its velocity is completely determined 
by exactly 4 parameters, computed in terms of the continued fraction 
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Figure 19. Collision pattern seen from the surface of 
one obstacle. Here, e = 2r/v\. 



expansion of Vijv\ — in the case where < i>2 < Vi, to which the 
general case can be reduced by obvious symmetry arguments. 

Assume therefore < vi < v \ with a — vijv\ Q. Henceforth, we 
set e = 2r\J\ + a 2 and define 

N(a, e) = inf {n > | d n (a) < e} , 



k(a, e) = — 



d 



N(a,e)-1 



dN («,£)(«) 



The parameters defining the collision pattern are A,B,Q — as they 
appear on the previous figure — together with an extra parameter 
S G {±1}- Here is how they are computed in terms of the continued 
fraction expansion of a = Vijv\. 



1 _ <^v(a, e )-i(QQ , fc(a,£)d JV(ctie )(a) 

1 £ " l " £ ' 



A(v,r) = l- dN(a f {a \ 

B(v,r) 
Q(v,r) = eq N ( a)£ )(a) , 
E(y,r) = (-l) N(a ' £ K 

The extra-parameter S in the list above has the following geometrical 
meaning. It determines the relative position of the closest and next 
to closest obstacles seen from the particle leaving the surface of the 
obstacle at the origin in the direction v. 
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The case represented on the figure where the closest obstacle is on 
top of the strip consisting of the longest particle path corresponds with 
S = +1, the case where that obstacle is at the bottom of this same 
strip corresponds with S = — 1. 

The figure above showing one example of collision pattern involves 
still another parameter, denoted Q' on that figure. 

This parameter Q' is not independent from A, B, Q, since one must 
have 

AQ + BQ' + A- B)(Q + Q') = 1 

each term in this sum corresponding to the surface of one of the three 
strips in the 3-term partition of the 2-torus. (Remember that the length 
of the longest orbit in the Blank-Krikorian theorem is the sum of the 
two other lengths.) Therefore 

1 — A(v, r) 

Once the structure of collision patterns elucidated with the help of 
the Blank-Krikorian variant of the 3-length theorem, we return to our 
original problem, namely that of computing the transfer map. 

In the next proposition, we shall see that the transfer map in a given, 
irrational direction « 6 S 1 can be expressed explicitly in terms of the 
parameters A, B, Q, E defining the collision pattern correponding with 
this direction. 

Denote 

K :=]0,l[ 3 x{±l} 

and let (A, B, Q, S) G K be the parameters defining the collision pat- 
tern seen by a particle leaving the surface of one obstacle in the direc- 
tion v. Set 

T A ,B, Q M h ') = (Q,h'-2E(1-A)) 

if 1 - 2A < Eh' < 1 , 
T A , B , Q ,v(h') = (Q',h' + 2E(1-B)) 

if -1 < Eh' < -1 + 25, 
T A B, Q Mh')= (Q' + Q,h' + 2E(A-B)) 

if -1 + 2B < Eh' < 1 - 2A . 

With this notation, the transfer map is essentially given by the explicit 
formula Ta,b,q,h, except for an error of the order 0(r 2 ) on the free-path 
length from obstacle to obstacle. 

Proposition 6.3 (Caglioti-Golse [TUl [H]) - One has 

T r (h',v) = T {AtB>Qi v )M (ti) + (O(r 2 ),0) 
in the limit as r — > + . 
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In fact, the proof of this proposition can be read on the figure above 
that represents a generic collision pattern. The first component in the 
explicit formula 



represents exactly rx the distance between the vertical segments that 
are the projections of the diameters of the 4 obstacles on the vertical 
ordinate axis. Obviously, the free-path length from obstacle to obstacle 
is the distance between the corresponding vertical segments, minus a 
quantity of the order 0(r) that is the distance from the surface of the 
obstacle to the corresponding vertical segment. 

On the other hand, the second component in the same explicit for- 
mula is exact, as it relates impact parameters, which are precisely the 
intersections of the infinite line that contains the particle path with the 
vertical segments corresponding with the two obstacles joined by this 
particle path. 

If we summarize what we have done so far, we see that we have 
solved our first problem stated at the beginning of the present section, 
namely that of finding a convenient way of coding the billiard flow in 
the periodic case and for space dimension 2, for a.e. given direction v. 

7. An ergodic theorem for collision patterns 

It remains to solve the second problem, namely, to find a conve- 
nient way of averaging the computation above so as to get rid of the 
dependence on the direction v. 

Before going further in this direction, we need to recall some known 
facts about the ergodic theory of continued fractions. 

The Gauss map 

Consider the Gauss map, which is defined on all irrational numbers 
in (0, 1) as follows: 



This Gauss map has the following invariant probability measure - 
found by Gauss himself: 



Moreover, the Gauss map T is ergodic for the invariant measure 
dg(x). By Birkhoff's theorem, for each / G L 1 (0, 1; dg) 




T : (0,l)\Q9ih4Ti 



X 



1 



11} e(0,l)\Q. 




1 



N-l „i 




N 



E/( T ^)^ / /(*)<&(*) 



a.e. in x G (0, 1) 
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as N — > +00. 

How the Gauss map is related to continued fractions is explained as 
follows: for 

a = [0; a , a u a 2 , • • •] = l — e (0, 1) \ Q 

a H ■ 

01 + . . . 

the terms afc(a) of the continued fraction expansion of a can be com- 
puted from the iterates of the Gauss map acting on a: specifically 

a k (a) = 

As a consequence, the Gauss map corresponds with the shift to the 
left on infinite sequences of positive integers arising in the continued 
fraction expansion of irrationals in (0, 1). In other words, 

T[0; a , ai, a 2 , • • •] = [0; a 1: a 2 , a 3 . . .] , 

equivalently recast as 

a n (Ta) = a n+ i(a) , n>0. 

Thus, the terms at(a) of the continued fraction expansion of any 
a G (0, 1) \ Q are easily expressed in terms of the sequence of iter- 
ates (T k a)k>o °f the Gauss map acting on a. The error d n (a) is also 
expressed in term of that same sequence (T k a)k>o, by equally simple 
formulas. 

Starting from the induction relation on the error terms 

d n +i{oi) = —a n (a)d n (a) + d n -i(a) , do(a) = 1 , d\(a) = a 
and the explicit formula relating a n (Ta) to a n (a), we see that 

ad n (Ta) = d n+ i(a) , n > . 
This entails the formula 

d n ( a ) = Y[T k a, n>0. 

k=0 

Observe that, for each 9 e [0, 1] \ Q, one has 

0-T0<\, 

so that 

d n (a) < , n>0, 



k > 
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which establishes the exponential decay mentionned above. (As a mat- 
ter of fact, exponential convergence is the slowest possible for the con- 
tinued fraction algorithm, as it corresponds with the rational approx- 
imation of algebraic numbers of degree 2, which are the hardest to 
approximate by rational numbers.) 

Unfortunately, the dependence of q n (a) in a is more complicated. 
Yet one can find a way around this, with the following observation. 
Starting from the relation 



we see that 



q n+1 (a)d n (a) + q n (a)d n+1 (a) = 1 , 

• d n (a)d n _i(at) 



g n (a)rf n _i(a) = ^(-1) 



n-j _ 



dj (a)dj-i(a) 

n—l 



= ^(-l) n ~ j X[T k - l aT k a. 

j=l k=j 

Usine; once more the inequality 9 ■ T9 < \ for 9 G [0, 1] \ Q, one can 
truncate the summation above at the cost of some exponentially small 
error term. Specifically, one finds that 



q n {a)d n -i{a) - ^ (-1] 



n-j 



:d n (a)d n -i(a) 



j=n-l 



dj(a)dj-i(a) 



n-l 



< T 



= q n (a)d n ^(a)- ^ (-1)^ J] T^aT'o 

j=n-l k=j 

More information on the ergodic theory of continued fractions can be 
found in the classical monograph [28J on continued fractions, and in 
Sinai's book on ergodic theory [30] , 

An ergodic theorem 

We have seen in the previous section that the transfer map satisfies 

T r (h',v) = T (AB ^ )M {h') + (O(r 2 ),0) as r - 0+ 

for each v G S 1 such that v 2 /vi G (0, 1) \ Q. 

Obviously, the parameters (A, are extremely sensitive to 

variations in v and r as r — > + , so that even the explicit formula for 
Ta,b,q,y,, is not too useful in itself. 

Each time one must handle a strongly oscillating quantity such as 
the free path length r r (x,v) or the transfer map T r (h',v), it is usually 
a good idea to consider the distribution of that quantity under some 
natural probability measure than the quantity itself. Following this 
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principle, we are led to consider the family of probability measures in 
(s,h) G R+ x [-1,1] 

S((s, h) - T r (h',v)) , 

or equivalently 

S((s,h) -T { A,B,Q,X)(v,r)(ti)) . 

A first obvious idea would be to average out the dependence in v of 
this family of measures: as we shall see later, this is not an easy task. 

A somewhat less obvious idea is to average over obstacle radius. 
Perhaps surprisingly, this is easier than averaging over the direction v. 

That averaging over obstacle radius is a natural operation in this 
context can be explained by the following observation. We recall that 
the sequence of errors d n (a) in the continued fraction expansion of an 
irrational a G (0, 1) satisfies 

ad n (Ta) = d n+ i(a) , n > , 

so that 

N(a, e) = inf {n > 1 1 d n (a) < e} 

is transformed by the Gauss map as follows: 

N(a,e) = N{Ta,e/a) + 1. 

In other words, the transfer map for the 2-dimensional periodic 
Lorentz gas in the billiard table Z r (meaning, with circular obstacles 
of radius r centered at the vertices of the lattice Z 2 ) in the direction v 
corresponding with the slope a is essentially the same as for the billiard 
table Z r / a but in the direction corresponding with the slope Tot. Since 
the problem is invariant under the transformation 

a i — > Ta , r i— > r/a 

this suggests the idea of averaging with respect to the scale invariant 
measure in the variable r, i.e. — on Rl. 

" r + 

The key result in this direction is the following ergodic lemma for 
functions that depend on finitely many d n s. 

Lemma 7.1 (Caglioti-Golse 01221 E]). For a G (0, 1) \ Q, set 

N(a, e) = inf {n > | d n (a) < e} . 
For each m > and each f G C(R™ +1 ) ; one has 

1 f 1/4 , ( ' d N ( at£ )(a) d/v(a, e )-m(aO\ de 

J 5 ■ ■ ■ ) > L m[J) 



| In 77I J v \ e e J e 

a.e. in a G (0,1) as r\ — > + ; where the limit L m (f) is independent of 



a. 
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With this lemma, we can average over obstacle radius any function 
that depends on collision patterns, i.e. any function of the parameters 
A,B,Q,E. 

Proposition 7.2 (Caglioti-Golse [II]). Let K = [0, l] 3 x {±1}. For 

each F G C(K), there exists C(F) G R independent of v such that 

1 C 1/2 dr 
— — / F(A(v,r),B(v,r),Q(v,r),Z(v,r))- ^ C(F) 

M^-fV) Jr, r 

for a.e. v G S 1 such that < V2 < v i in the limit as rj — > + . 

Sketch of the proof. First eliminate the £ dependence by decomposing 

F(A, B, Q, S) = F+(A, B, Q) + 5, Q) . 

Hence it suffices to consider the case where F = F(A, B,Q). 
Setting a — Vijv\ and e = 2r/vi, we recall that 

At \ ■ f f d N(a,e){ a ) 

A{v,r) is a function of , 

B{v,r) is a function of — - — and — - — . 

As for the dependence of F on Q, proceed as follows: in F(A, B, Q), 
replace Q(v,r) with 

g / ^ N{a,e)-j d N(a,s) («)^V(q, £ )-l («) 

v ' j=N(a,e) — l 

at the expense of an error term of the order 

0(modulus of continuity of F{2~ m )) — >• as / —> , 

uniformly as £ — » + . 

This substitution leads to an integrand of the form 

/ djv(a, £ )(a) C? A r( ai£ )_ m _ 1 (Q;) 
J I ) • • • ) 

\ £ e 

to which we apply the ergodic lemma above: its Cesaro mean converges, 
in the small radius limit, to some limit C m (F) independent of a. 
By uniform continuity of F, one finds that 

\C m (F) - C m >{F)\ = O (modulus of continuity of F{2- mVm ')) 

(with the notation m V m! = max(t>, v')), so that C m (F) is a Cauchy 
sequence as m — ► oo. Hence 

C m (F) — > C{F) as m — ► oo 
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and with the error estimate above for the integrand, one finds that 



With the ergodic theorem above, and the explicit approximation of 
the transfer map expressed in terms of the parameters (A, B, Q, S) 
that determine collision patterns in any given direction v, we easily 
arrive at the following notion of a "probability of transition" for a 
particle leaving the surface of an obstacle with an impact parameter h! 
to hit the next obstacle on its trajectory at time s/r with an impact 
parameter h. 



Theorem 7.3 (Caglioti-Golse, [101111]). For each h' E [-1,1], there 
exists a probability density P(s,h\h') on R + x [—1,1] such that, for 
each f G C(R+ x [-1, 1]), 



a.e. in v G S 1 as r] — > + . 

In other words, the transfer map converges in distribution and in the 
sense of Cesaro, in the small radius limit, to a transition probability 
P(s, h\h') that is independent of v. 

We are therefore left with the following problems: 

a) to compute the transition probability P(s,h\h') explicitly and dis- 
cuss its properties, and 

b) to explain the role of this transition probability in the Boltzmann- 
Grad limit of the periodic Lorentz gas dynamics. 



Most unfortunately, our argument leading to the existence of the 
limit C(F), the core result of the previous section, cannot be used for 
computing explicitly the value C(F). Indeed, the convergence proof is 
based on the ergodic lemma in the last section, coupled to a sequence 
of approximations of the parameter Q in collision patterns that involve 
only finitely many error terms d n (a) in the continued fraction expansion 
of a. The existence of the limit is obtained through Cauchy's criterion, 
precisely because of the difficulty in finding an explicit expression for 
the limit. 




as rj 



0+ 



□ 




8. Explicit computation 

OF THE TRANSITION PROBABILITY P(s,h\h' 
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Nevertheless, we have arrived at the following expression for the 
transition probability P(s,h\h'): 

Theorem 8.1 (Caglioti-Golse [M1III]). The transition probability den- 
sity P(s, h\h') is expressed in terms of a = — h'\ and b = + h'\ 
by the explicit formula 

P(s,h\ti) = -|_[(( s -I sa )A(l + iso)-(lV(|s + |s6)) 



n 2 sa 



+ ((s- isa)Al-((|s + |sfe)v(l-|sa)) + 

+saA|l — s|l s< i + (sa — 11 — s|) + , 

with the notations x Ay = min(a;, y) and x V y = max(i, y). 
Moreover, the function 

(s, fr') i — > (l -\- s)P(s, h\ti) belongs to L 2 (R + x [-1, l] 2 ) . 

In fact, the key result in the proof of this theorem is the asymptotic 
distribution of 3-obstacle collision patterns — i.e. the computation of 
the limit C(f), whose existence has been proved in the last section's 
proposition. 

Theorem 8.2 (Caglioti-Golse plj). Define K = [0, l] 3 x {±1}; then, 
for each F eC(K) 

1 r 1/4 dr 
—— / F{(A,B,Q,Ti)(v,r))-^C(F) 
\\nr]\ J v r 

= / F{A, B, Q, T)dm{A, B, Q, E) a.e. m v G S 1 
Jk 

as r] — > + , where 

dm(A, B, Q, E) = dm (A, B, Q) g> |(5 s =i + fe=-i) , 

dm (A, B,Q) = ^|1 0<j 4 <1 1 <b<i_a1o<q< 



l2 dAdBdQ 

2-A-B I - A 



Before giving an idea of the proof of the theorem above on the dis- 
tribution of 3-obstacle collision patterns, it is perhaps worthwhile ex- 
plaining why the measure m above is somehow natural in the present 
context. 

To begin with, the constraints < A < 1 and < B < 1 — A have 
an obvious geometric meaning (see figure 18 on collision patterns.) 
More precisely, the widths of the three strips in the 3-term partition 
of the 2-torus minus the slit constructed in the penultimate section (as 
a consequence of the Blank-Krikorian 3-length theorem) add up to 1. 
Since A is the width of the strip consisting of the shortest orbits in 
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the Blank-Krikorian theorem, and B that of the strip consisting of the 
next to shortest orbits, one has 

< A + B < 1 

with equality only in the exceptional case where the orbits have at 
most 2 different lengths, which occurs for a set of measure in v or r. 
Therefore, one has 

< B(v, r) < 1 - A(v, r) , for a.e. r e (0, \) . 

Likewise, the total area of the 2-torus is the sum of the areas of the 
strips consisting of all orbits with the 3 possible lengths: 

l = QA + Q'B + (Q + Q')(l - A - B) = Q(l - B) + Q'(l - A) 

>Q(2-A-B) 

as Q' > Q (see again the figure above on collision patterns.) 
Therefore, the volume element 

dAdBdQ 
I- A 

in the expression of dm imples that the parameters A, - - or 
equivalently B mod. 1 — A — and Q are independent and uniformly 
distributed in the largest subdomain of [0, l] 3 that is compatible with 
the geometric constraints. 

The first theorem is a consequence of the second: indeed, P(s, h\h')dsdh 
is the image measure of dm(A, B, Q, E) under the map 

K 9 (A, 5, Q, £) ^ T (ABiQiS) (/i>) . 

That (1 + s)P(s, h\h!) is square integrable is proved by inspection - 
by using the explicit formula for P(s, h\h'). 

Therefore, it remains to prove the second theorem. 
We are first going to show that the family of averages over velocities 
satisfy 

F{A{v,r),B{v,r),Q(v,r),Y.{v,r))dv 

v\=l 

-f / F(A,B,Q,E)dm(A,B,Q,E) 
Jk 

as r -> 0+ for each F E C b (K). 

On the other hand, because of the proposition in the previous section 

1 r 1/2 dr 

J F(A(v, r), B(v, r), Q(v, r), E(v, r)) - -> C(F) 

for a.e. v G S 1 such that < v 2 < v i in the limit as rj — > + . 
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Since we know that the limit C(F) is independent of v, comparing 
the two convergence statements above shows that 

C{F)= [ F(A, B, Q, E)dm(A, B, Q, S) . 

Therefore, we are left with the task of computing 

lim / F(A(v, r), B(v, r), Q(v , r), E(v, r))dv . 

r^0+ J M=l 

The method for computing this type of expression is based on 

a) Farey fractions (sometimes called "slow continued fractions"), and 

b) estimates for Kloosterman's sums, due to Boca-Zaharescu [3]. 
To begin with, we need to recall a few basic facts about Farey frac- 
tions. 

Farey fractions 

Put a filtration on the set of rationals in [0, 1] as follows 

Fq = {? I o < v < q < Q , g.c.d.(p, q) = 1} 
indexed in increasing order: 

Pi 1 

= - < 7i < • • • < H = < ■ ■ ■ < 1 V {Q) = t = 1 

where if denotes Euler's totient function: 



p prime 
p|n 

An important operation in the construction of Farey fractions is the 
notion of "mediant" of two fractions. Given two rationals 

7 = - and 7 = - 
q q 

with 

0<p<q, 0<p<q, and g.c.d.(p, q) = g.c.d.(p, q) = 1 , 
their mediant is defined as 

mediant = 707:= G (7, 7) . 

q + q 

Hence, if7 = 2<^ = | adjacent in jFg, then 
aq — aq = 1 and q + q > Q . 

Conversely, q, q are denominators of adjacent fractions in Tq if and 
only if 

< g, q<Q, q + q>Q, g.c.d.(q, q') = 1 . 
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143525341 

\ /l\ ! /!\ ! /II ! /!\ ! 

Figure 20. The Stern-Brocot tree. Each fraction 7 on 
the ra-th line is the mediant of the two fractions closest 
to 7 on the 71 — 1-st line. The first line consists of and 
1 written as = ° and 1 = \. Each rational in [0, 1] is 
obtained in this way. 



Given a G (0, 1)\Q and Q > 1, there exists a unique pair of adjacent 
Farey fractions in JF Q , henceforth denoted 7(0;, Q) and j(a, Q), such 
that 

At this point, we recall the relation between Farey and continued 
fractions. 

Pick < e < 1; we recall that, for each a G (0, 1) \ Q, 
N(a, e) = min{n G N | d n (a) < e} , d n (a) = dist(q n (a)a, Z) . 
Set Q = [1/e], and let 

g(a, Q) Q) 

with g.c.d.(p(a, Q)), q(a, Q)) = g.c.d.(p(a, Q)), q(a, Q)) = 1 be the two 
adjacent Farey fractions in Tq surrounding a. Then 

a) one of the two integers q(a, Q) and q(a, Q) is the denominator of 
the N(a, e)-th convergent in the continued fraction expansion of a, i.e. 

QN( a ,e)(ot), and 
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b) the other is of the form 

mq N ( a ,e) + qN( a ,e)-i , with < m < a N ^ £) {a) , 
where we recall that 

a = [0;a 1 ,a 2 ,...\ = 



1 

a + 



ax + . . . 

Setting a = V2/V1 and e = 2r/vi, we recall that, by definition 
Q(v,r) = eq N(a>£) (a) G {eq(a, Q),eq(a, Q)} with Q = [l/e\ , 
and we further define 

D(v,r) = dtf (a>e )/e = dist(J<3(i>, r)a, Z)/e , 

and 

Q(i>,r) = Q) if 5jv(a,e)(a) = ?(<*, Q) , 
Q(v,r) = eq(a, Q) if q N ^ £ )(a) = q(a, Q) . 
Now, we recall that A(v,r) = 1 — D(v,r); moreover, we see that 



B(v, r) 



1 - d N{a)£) -i(a)/s 



D(v, r) 

1 - djv(a,e)-i(a)/£ mod. £>(w,r) 



D(f , r) 



= 1 — &\st(\Q(v, r)ot, Z)/e mod. r) . 
To summarize, we have 

F(A(v,r),B(v,r),Q(v,r)) = G(Q(v,r),Q(v,r), D(v,r)) 
and we are left with the task of computing 

lim / G(Q(v,r),Q(v,r), D(v,r))dv 
r ^o+ J s i + 

where is the first octant in the unit circle. The other octants in the 
unit circle give the same contribution by obvious symmetry arguments. 
More specifically: 

Lemma 8.3. Let a G (0, 1) \Q 7 and let | < a < | be the two adjacent 
Farey fractions in Tq surrounding a, with Q = [l/e]. Then 
a) if v - <a<^, then 

/ J q — q > 

Q(v,r)=eq, Q(v,r)=eq, D(v,r) = Maq - p) ; 
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b) ifVf <a< § ; then 

Q(v,r)=eq, Q(v,r) = eq, D(v , r) = i(p - acq) ; 

c) z/£±^ < a < tfien 

/ J q — q ' 

Q(v,r) = e(qAq) , Q(v,r) = s(qVq) , D(v, r) = dist(±Q(v, r)a, Z) . 
Therefore, assuming for simplicity that 

y, z) = o(:r, y)H'(z) and e = 1/Q , 

one has 



E r ^ 9{%^)H'{Q{qa-v))da 

n< O <? n4-n 



G{Q{v,r),Q{v,r),D(v,r))dv 

*(p-e)/« 

0<9,q<S< 9 + tj 
g.c.(i.(gj)=l 

+ three other similar terms 



0<q,q<Q<q + q 
g.c.d.(q,q) = l 



+ three other similar terms. 
Then, everything reduces to computing 



^ 0<n,<5<S<<7 + 9 \ ^ 



3<<j,9<S<g+<? 

g.c.d.(q,q) = l 



We conclude with the following 

Lemma 8.4 (Boca-Zaharescu [3]). For ip e C C (R 2 ) ; one /ias 



Q2 £ ^(q'Q/ ^ 



* 0< 9 ,9<Q< 9 + <3 ^ ^0<I,J/<1<I+J 

g.c.d.(g,q)=l 

in the limit as Q — > oo. 



^(x, y)dxdy 



This is precisely the path followed by F. Boca and A. Zaharescu to 
compute the limiting distribution of free path lengths in [3] (see Theo- 
rem 4.6); as explained above, their analysis can be greatly generalized 



in order to compute the transition probability that is the limit of the 
transfer map as the obstacle radius r — > + . 
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9. A KINETIC THEORY IN EXTENDED PHASE-SPACE 
FOR THE BOLTZMANN-GRAD LIMIT 
OF THE PERIODIC LORENTZ GAS 

We are now ready to propose an equation for the Boltzmann-Grad 
limit of the periodic Lorentz gas in space dimension 2. For each r e 
]0, |[, denote 

B r : 9 (x, v) i— > B r (x, v ) = (x + T r (x, v)v, 1Z[x + T r (x, v)v]v) E 
the billiard map. For (x ,v ) G T+, set 

(x n ,v n ) = B™(x ,v ) 

and define 

b?(x,v) = (A,B,Q,Ti)(v n ,r), neN*. 
Henceforth, for each n > 1, we denote 

/C n := R 2 x S 1 x R + x [-1, 1] x K n . 

We make the following asymptotic independence hypothesis: there 
exists a probability measure II on R + x [— 1, 1] such that, for each n > 1 
and each \1/ e C(A^ n ) with compact support 

lim / <&(x,v,rT r {f,v),h r (f,v 1 ),b l r ,...,b^)dxdv 



ZrXS 1 

*$?(x, v, t, h, , /3 n )dxdvdH(T, h)dm(Pi) . . . dm(f3 n 



Qn 



where 

0o,t>o) = (x- r r (x, -v)v,v) , and /i r (^,wi) = sm(n xl ,Vi) 

and m is the probability measure on K obtained in Theorem 8.2| 
If this holds, the iterates of the transfer map T r are described by the 

Markov chain with transition probability P(s, h\h'). 

This leads to a kinetic equation on an extended phase space for the 

Boltzmann-Grad limit of the periodic Lorentz gas in space dimension 

2: 

F{t, x, v, s, h) = 

density of particles with velocity v and position x at time t 
that will hit an obstacle after time s, with impact parameter h. 

Theorem 9.1 (Caglioti-Golse [TO] HI]). Assume (H), and let f in > 
belong to C C (R 2 x S 1 ). Then one has 

/ F(", -, •, s, h)dsdh in L°°(R + x R 2 x S 1 ) weak-* 
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in the limit as r — > + 7 where F = F(t,x,v, s, h) is the solution of 

(d t +v ■ W x - d s )F(t, x, v, s, h) 
i 

P(s, h\h')F(t, x, R[n - 2 arcsin(/i')K 0, , 

i 

y P(T,h\h')dh'dr, 

with (x,v, s, /i) running through R 2 x S 1 x R^x] — 1, 1[. T/ie notation 
R[9] designates the rotation of an angle 9. 

Let us briefly sketch the computation leading to the kinetic equation 
above in the extended phase space Z = R 2 x S 1 x R + x [—1,1]. 

In the limit as r — > + , the sequence (b™(x, v)) n >i converges to a 
sequence of i.i.d. random variables with values in K = [0,1] x {±1}, 
according to assumption (H). 

Then, for each s > and h G [—1,1], we construct a Markov chain 
(•Sri, h n ) n >i with values in R + x [—1, 1] in the following manner: 

(s n ,h n ) = T bn (h n _ 1 ) , n>l. 

Now we define the jump process (X t , V t , S t , H t ) starting from (x, v, s, h) 
in the following manner. 

First pick a trajectory of the sequence b = (b n ) n >i, then, for each 
s > and each h G [—1, 1], set 

(s , ho) = (s,h) . 

Define then inductively s n and h n for n > 1 by the formula above, 
together with 

0"n = so + • • • + s n _i , n > 1 , 

and 

v n = R[2 arcsin(/i n _i) — n]v n -i , n > 1 . 

With the sequence (v n , s n , h n ) n >i so defined, we next introduce the 
formulas for (X t , V t , S t ,H t ): 

• While < t < t, we set 

X t (x, v, s, h) — x + too , S t (x,v,s,h) —s — t, 

V t (x, v, s, h) — v , H t (x,v,s,h) —h. 
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• For a n < t < cr n+ i, we set 

X t (x, v,s,h) = x + (t - a n )v n , 

V t (x,v, s, h) = v n , 

T t (x,v,s,h) = <7 n+ i - i, 
H t (x,v, s, h) = h n . 

To summarize, the prescription above defines, for each t > 0, a map 
denoted T t : 



2xK N 9 (x, u, s, fc, b) ^ T t (x, r, h) = (X t , V t , S t , H t ) e Z 



that is piecewise continuous in t G R + . 

Denote by f m = f m (x,v, s,h) the initial distribution function in 
the extended phase space Z, and by x = x{ x i v i s > ^) an observable - 
without loss of generality, we assume that x £ C^°(Z). 

Define f(t, ■,-,-,■) by the formula 



where E designates the expectation on trajectories of the sequence of 
i.i.d. random variables b = {b n ) n >\. 

In other words, f(t, •, •, •, •) is the image under the map T t of the 
measure Prob(cib)/ m (a;, u, r, h), where 



g(t, x, v, s, h) = E[l t<s x(T t (x, v, s, h)))] + E[l s<t x(T t (x, v, s, h)))} . 



If s > t, there is no collision in the time interval [0, t] for the trajectory 
considered, meaning that 




Prob(db) = Yl dm{b n ) . 



n>l 

Set g(t, x, v, s, h) = E[x(T t (x, v, s, h)))]; one has 



T t (x, v, s, h) — (x + tv, v, s — t,h) . 



Hence 



E[l t <sX( T t(x, v, s, h)))] = x(x + tv, v,s-t, h)l t<s . 



On the other hand 



E[l s<tX (T t (x, v, s, h))} = E[l s<tX (T {t _ s y T s+0 (x, v, s, h))} 
= E[l s<tX (T {t _ s) _ (x + sv, K[A(h)]v, si, h 1 ))} 

with (si, hi) = T bl (h) and A(h) = 2arcsin(/i) — it. 
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Conditioning with respect to (si, hi) shows that 
E[l s<tX {T t {x,v,s,h)))} 

= E[l s<t E[x(T (t _ s) _ (x + sv, K[A(h)]v, s u hi))\s u h}] , 

and 

E[x(T ( t_ s) _ (x + sv,K[A(h))v,Si,hi))\si,hi) 
= g(t — s, x + sv, lZ[A{h)}v, Si, hi) . 

Then 

E[l s<t E[ X (T (t _ s) ^(x + sv, TZ[A(h)]v, si, hi))\s u h)] 
= l s<t J g(t-s,x + sv,K[A(h))v,T bl {h)))dm{bi) 

— l s<t J g(t — s,x + sv,1Z[A(h)}v,si,hi))P(si,hi\h)dsidhi . 
Finally 

g(t, x, v, s, h) = x{x + tv, v,s-t, h)l t<s 

+ l s<i J g(t — s,x + sv,lZ[A(h))v, Si, h 1 )]P(s 1 , hi\h)dsidhi . 

This formula represents the solution of the problem 

(d t - v ■ + d s )g = , t,s>0, x e R 2 , s e S 1 , \h\ < 1 



g(t, x, s, 0, h) = // P(si, hi\h)g(t, x,v, Si, hi)dsidhi , 

J 7r* x]-i,i[ 
g\ t=0 = x- 

The boundary condition for s = can be replaced with a source term 
that is proportional to the Dirac measure S s= q: 

{d t - v ■ V x + d s )g = 5 S=0 // P{si,h 1 \h)g{t,x,v,si,hi)ds l dh l , 

J 7r* x]-i,i[ 

One concludes by observing that this problem is precisely the adjoint 
of the Cauchy problem in the theorem. 

Let us conclude this section with a few bibliographical remarks. 

Although the Boltzmann-Grad limit of the periodic Lorentz gas is a 
fairly natural problem, it remained open for quite a long time after the 
pioneering work of G. Gallavotti on the case of a Poisson distribution 
of obstacles [TBI US]- 

Perhaps the main conceptual difficulty was to realize that this limit 
must involve a phase-space other than the usual phase-space of kinetic 



PERIODIC LORENTZ GAS 



59 



theory, i.e. the set R 2 x S 1 of particle positions and velocities, and to 
find the appropriate extended phase-space where the Boltzmann-Grad 
limit of the periodic Lorentz gas can be described by an autonomous 
equation. 

Already Theorem 5.1 in [9] suggested that, even in the simplest prob- 
lem where the obstacles are absorbing — i.e. holes where particles 
disappear forever, — the limit of the particle number density in the 
Boltzmann-Grad scaling cannot described by an autonomous equation 
in the usual phase space R 2 x S 1 . 

The extended phase space R 2 x S 1 x R + x [—1,1] and the structure 
of the limit equation were proposed for the first time by E. Caglioti 
and the author in 2006, and presented in several conferences — see 
for instance [23] ; the first computation of the transition probability 



P(s,h\h') (Theorem 8.1), together with the limit equation (Theorem 



9.1) appeared in (TU] for the first time. However, the theorem concern- 



ing the limit equation in [TU] remained incomplete, as it was based on 
the independence assumption (H). 

Shortly after that, J. Marklof and A. Strombergsson proposed a com- 
plete derivation of the limit equation of Theorem 9A in a recent preprint 
|32j . Their analysis, establish the validity of this equation in any space 
dimension, using in particular the existence of a transition probabil- 
ity as in Theorem 8A_ in any space dimension, a result that they had 
proved in an earlier paper [31] . The method of proof in this article [31] 
avoided using continued or Farey fractions, and was based on group ac- 
tions on lattices in the Euclidian space, and on an important theorem 
by M. Ratner implying some equidistribution results in homogeneous 
space. However, explicit computations (as in Theorem 8.1 of the tran- 
sition probability in space dimension higher than 2 seem beyond reach 
at the time of this writing — see however [33J for computations of the 
2-dimensional transition probability for more general interactions than 
hard sphere collisions. 

Finally, the limit equation obtained in Theorem 9.1 is interesting in 
itself; some qualitative properties of this equation are discussed in [TT] . 



Conclusion 

Classical kinetic theory (Boltzmann theory for elastic, hard sphere 
collisions) is based on two fundamental principles 

a) deflections in velocity at each collision are mutually independent 
and identically distributed 

b) time intervals between collisions are mutually independent, inde- 
pendent of velocities, and exponentially distributed. 
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The Boltzmann-Grad limit of the periodic Lorentz gas provides an 
example of a non classical kinetic theory where 

a') velocity deflections at each collision jointly form a Markov chain; 

b') the time intervals between collisions are not independent of the 
velocity deflections. 

In both cases, collisions are purely local and instantaneous events: 
indeed the Boltzmann-Grad scaling is such that the particle radius 
is negligeable in the limit. The difference between these two cases is 
caused by the degree of correlation between obstacles, which is maxi- 
mal in the second case since the obstacles are centered at the vertices of 
a lattice in the Euclidian space, wheras obstacles are assumed to be in- 
dependent in the first case. It could be interesting to explore situations 
that are somehow intermediate between these two extreme cases — for 
instance, situations where long range correlations become negligeable. 

Otherwise, there remain several outstanding open problems related 
to the periodic Lorentz gas, such as 

i) obtaining explicit expressions of the transition probability whose 
existence is proved by J. Marklof and A. Strombergsson in |31j . in all 
space dimensions, or 

ii) treating the case where particles are accelerated by an external 
force — for instance the case of a constant magnetic field, so that the 
kinetic energy of particles remains constant. 
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